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ABSTRACT 

The average capture cross section for gold ct c at a neutron energy of 24 keV has 
been obtained by analysis of published sphere transmission experiments that use Sb-Be 
source neutrons. The analysis compares results using merged s-wave spin -state 
statistical data and separated spin-state statistical data. Values of a c of 725±30 mb 
are obtained using separated spin-state data and 660±30 mb using merged spin-state 
data. The value of 725±30 mb is in agreement with the larger of recent gold activation 
cross section measurements at 24 keV, which have a range of from 640±25 mb to 
713±40 mb for the average capture cross section 


ANALYSIS OF GOLD SPHERE TRANSMISSION EXPERIMENTS AT 24-keV 
NEUTRON ENERGY USING SPIN-DEPENDENT s-WAVE STATISTICS 

by Thor T. Semler 
Lewis Research Center 

SUMMARY 

A Monte Carlo analysis of spherical shell transmission experiments for gold at a 
neutron energy of 24 keV is presented. In the analysis, published spin -dependent neutron 
strength functions are used to generate s-wave cross sections. 

A total of 10 cross sections are generated. A choice of the best sample of 10 4 cross 
sections, as determined by the Central Limit Theorem is made. These cross sections 
are used as input to the transmission program. 

The potential cross section a is chosen as a parameter to conserve the independ- 
ently measured value of the total cross section ^ for gold of 13. 67 barns. For a given 
value of cr j., several values of the average p-wave capture cross section a are 

p c p 

introduced, and values of the neutron transmission are computed for each experimental 

shell. The value of a compatible with the experimentally measured transmission is 

P 

found by interpolation of several Monte Carlo solutions. This procedure is repeated for 
several values of a ^ The final value of a~ is obtained as the value compatible with 
P P 

the measured transmission for a particular shell, and the total cross section since 

°tot = °pot + a s + a c + o* c (where cT“ and are, respectively, the average s-wave 

s s p s s 

scattering and capture cross sections). 

Thus one obtains from the Monte Carlo analysis and the shell transmission experi- 
ments a set of cross sections which match the shell transmission and conserve the value 
of the total cross section. 

A value of o~^ at 24 keV for gold of 725±30 millibarns is obtained using separated 
s-wave spin-state data compared with a value of 660±30 millibarns using merged s-wave 
spin-state data. This value is in agreement with the larger values of recent gold acti- 
vation cross section measurements at 24 keV, which range from 640±25 to 713±40 milli- 
barns for the average capture cross section 
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INTRODUCTION 

Neutron capture cross section measurements in the keV energy region are of con- 
siderable interest to fields as diverse as cosmology and reactor design. Of particular 
interest are precise neutron capture cross sections for gold which could serve as a con- 
venient secondary standard in the measurement of other isotopic capture cross sections. 
The status of the gold capture cross sections has been discussed recently by Cox (ref. 1) 
and by Poenitz (ref. 2). Attempts have been made to explain the discrepancies and to 
reconcile the various measurements in the keV region by Bogart (ref. 3) and by Poenitz, 
Kompe, and Menlove (ref. 4) who have discussed the persisting difficulties. 

Most of the measurements below 200 keV rely on a normalization to activation cross 
sections at 30 keV or to sphere transmission results using antimony -beryllium (Sb-Be) 
source neutrons of nominally 24 keV. A summary of these cross sections for gold has 
been presented by Ryves et al. (ref. 5) who notes that recent activation cross sections 
amended to correspond to 24. 8 keV vary from 713±40 millibarns of Harris et al. (ref. 6) 
to 679±70 millibarns of Poenitz (ref. 2) to 640±25 millibarns of Ryves et al. (ref. 5). 

The shell transmission measurements provided values of 532±60 and 605±30 milli- 
barns (refs. 7 and 8, see table I). However, the corrections for resonance self- 
protection that had been applied to the sphere transmission experiments (ref. 7) were 
found to be inadequate. A Monte Carlo interpretation of these experiments (reported in 
refs. 9 and 10) using merged spin s-wave statistics to generate cross sections at 24 keV, 
provided larger values of cross sections of 635±50 millibarns (ref. 7) and 660±50 milli- 
barns (ref. 8). 


TABLE I. - PREVIOUS ANALYSES OF GOLD-197 
SPHERE TRANSMISSION EXPERIMENTS 


Experiment 

Data from - 


Original 

source 

Bogart and Semler 
(ref. 9) 

This work 


Capture cross sections, c 

, mb 
c’ 

a Au I 

532±60 

620±50 

710±31 

a Au n 

518±90 

635±60 

737±33 

b Long counter 

570±30 

640±60 

774±31 


a Ref. 7. 
b Ref . 8. 
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With the recent availability of separated spin -dependent s-wave statistics for gold 
(ref. 11), the Monte Carlo analysis of the sphere transmission experiments has been 
remade in an attempt to improve the precision of the average gold capture cross section 
at 24 keV. 


SPHERE TRANSMISSION EXPERIMENT AND ANALYSIS 

The experiments to be analyzed (refs. 7 and 8) used a monoenergetic Sb-Be neutron 
source inside of a spherical shell of the material whose cross section is to be measured. 
A flat-response neutron detector or water bath measures the number of neutrons emer- 
gent for the shell surrounding the source and the bare source. The observed values of 
transmission are corrected for finite source size and other effects, including the small 
group of higher energy neutrons emitted by the Sb-Be source. The corrected values of 
transmission are the ratios of the number of neutrons escaping the spherical shell of 
material to those emitted by a point source of 24 keV neutrons centrally located in the 
spherical shell of material. 

From these values of transmission, a value for the average capture cross section 
can be obtained by an analysis due to Bethe et al. (refs. 12 to 14) which assumes cross 
sections to vary smoothly with neutron energy. Although the assumption of smoothly 
varying cross sections may be valid in the MeV energy region, at lower energies the 
assumption is questionable because of the significant peaking of individual resonances. 
The difficulties encountered when analysis of references 12 to 14 is applied to transmis- 
sion experiments at 24 keV have been illustrated in reference 9. 

The Monte Carlo analyses made herein follow the multiple scatterings of the neu- 
trons in the material, and use directly the pointwise cross sections in the unresolved 
energy region generated by level statistics obtained from the resolved region. 


MONTE CARLO TRANSMISSION CALCULATION 

A Monte Carlo spherical transmission code is used to analyze the gold shell trans- 
mission experiments of Schmitt and Cook (ref. 7) and those of Belanova et al. (ref. 8). 
The Monte Carlo code assumes a point isotropic source of neutrons located at the center 
of a hollow spherical shell of material. The collision trajectories are followed until the 
neutron is absorbed or has been transmitted from the shell. On elastic collision, the 
neutron is scattered isotropically in the center -of -mass system, and for a heavy nucleus, 
such as gold, it is a reasonable simplification of the calculations if elastic scattering is 
taken to be isotropic in the laboratory system as well. 
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Start 



Figure L - Schematic description of Monte Carlo sphere transmission calculation. 


A schematic description of the Monte Carlo sphere transmission calculation is shown 
in figure 1 and is further explained as follows: 

(1) A cross section is chosen. (The generation of cross section is described in the 
next section. ) 

(2) A path length A r is chosen from an exponential distribution by the following 
equation: 

-In N r 
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where N r is a random number uniformly distributed on the open interval (0, 1) 
and S t is the macroscopic total cross section. 

(3) If A r terminates outside the shell, the neutron is tabulated as a transmission 
and a return is made to step 1. 

(4) If X r terminates inside the spherical shell, the direction cosines of an isotrop- 
ically emitted neutron are computed, and a new position of the neutron interaction is 
computed (ref. 15). 

The cross sections are selected, and, if the collision is an absorption, the history 
is terminated. If the collision is a scattering, a return is made to step 1. The final 
output includes the transmission and the average cross sections for scattering and absorp- 
tion. 


STOCHASTIC CROSS SECTION GENERATION 

The method of stochastic simulation is used to generate cross sections at 24 keV in 
order to accurately represent the behavior of the neutron trajectories in the spherical 
shell. This is done by applying the statistical distributions obtained from neutron reson- 
ance spectroscopy in the resolved energy region to neutron cross sections at 24 keV. 

1Q7 

Desjardins et al. (ref. 16), in their measurements of the Au resonance param- 
eters, were unable to assign J, the spin quantum number of the compound nucleus values 
to many of the resonances they resolved. Therefore only average values of the resonance 
parameters for the merged spin states were reported. The best -fit value of the merged 
s-wave strength function S Q was 


<gr°> 4 

= J 1 - = (1. 5±0. 3)xl0" 4 

° D 

where g is a statistical spin factor equal to (2J + 1)/ (41 + 2) and I is the spin quantum 

number of the nucleus investigated. The integral distribution of the reduced width neutron 

amplitudes r° showed good agreement with a Porter -Thomas distribution. The merged 

spin average level spacing D was found to be 15. 8 eV below 1000 eV, and 16. 8 eV was 

suggested for the entire measured range to 1600 eV. However, when the integral distri- 

197 

button of gold- 197 (Au ) level spacings was compared with the distributions, for a 
random spacing of levels and with a single-species population Wigner distribution; the 
agreement is poor (ref. 17). Thus there was no agreement of the level spacing results 
of Desjardins et al. and these theoretical curves. Since gold has a nuclear spin 1 = 3/2, 
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the J = 1 levels axe assumed to have a Wigner single-species population distribution, 
which is uncorrelated in position with respect to the J = 2 levels that are also assumed 
to be Wigner single -species population distribution. If the (2J + 1) spacing density 
applies, 3/8 of the levels should have J = 1 and the remaining 5/8 levels have J = 2. 
The Wigner two-species population curves, D x = 5D 2 /3, gave good agreement when com- 
pared with the initial distribution; in fact, if Dj = equally good fits can be obtained. 
The average parameters measured by Desjardins et al. of Columbia have been adjusted 
to 24 keV and are shown in table n . 

TABLE H. - GOLD -197 STATISTICAL PARAMETERS USED IN 


CROSS SECTION GENERATOR 


— 

Source 

Type of 
spins 

Spin 

quantum 
number , 

I 

Statistical 

spin 

factor, 

g 

y width, 

r r’ 

meV 

Neutron 

width, 

r n> 

meV 

Average 

level 

spacing, 

D, 

eV 

Ref. 16 

Merged 

1 and 2 

0. 500 

170 

780 

16. 8 

Ref. 11 

Merged 

1 and 2 

. 500 

125 

930 

15. 5 


Separated 

1 

.375 

125 

930 

44. 0 



2 

.625 

125 

930 

24. 0 


More recently, Julien et al. (ref. 11) of Saclay have measured the spin -dependent 
statistical parameters of gold. Their results were different than those of Desjardins 
et al. The Saclay results were more in agreement with a value of D^ = 2. 2 T>^. They 
found the reduced level width distribution functions as measured to agree with the theo- 
retical Porter-Thomas distribution. They also found far better agreement between their 
results and the Wigner two-species population distribution than with a Wigner one -species 
population distribution. The average separated spin statistical parameters, as well as 
the merged values generated from the Saclay data, have been adjusted to 24 keV and are 
shown in table II . 

With the statistical information obtained at low neutron energies, a series of "mock" 
or "pseudo” cross sections for the unresolved region are generated by the method of 
stochastic simulation. By making use of an adequate cross section set, as determined 
by the Central Limit Theorem, and a large energy spread for all resonances chosen 
(E q ±100 eV), a so-called resonance background is provided without the problem of the 
distribution of next nearest levels. The contribution of nearby resonances to the cross 
section is considered if the resonance is within 100 eV of the energy of the sampled cross 
section's energy. The use of a mock cross section set, decided on by these statistical 
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methods, ensures an average not subject to large fluctuations. Since both sets of s-wave 
capture cross sections are computed from independently chosen random number 
sequences, the uncorrelated aspect of the separate J -state resonances is assured. 

Since calculation indicates that the peak values of the Doppler -broadened p-wave 

197 

scattering resonances for Au at 24 keV are only of the order of a few tenths of a barn, 
the p-wave scattering is explicitly neglected, and any contribution to the total cross sec- 
tion by p-wave scattering is subsumed into the potential cross section. 

In as much as both the Columbia and Saclay (refs. 7 and 8) data indicate that the 

o 

distribution of the reduced scattering width amplitudes is Porter-Thomas, that is x , 
the values for r n have been chosen from a x distribution with one degree of freedom. 
The selection method is detailed in the appendix. Briefly, the values of T n (J = 1) and 
r n (J = 2) are computed for 24-keV neutrons. Next, the cumulative distribution function 
of r n (J = l)/r n (J = 1) is divided into 10 intervals of equally probable values. The 
average value of r n / T n is found in each interval, and the corresponding value of T n 
chosen. Thus, one is able to represent the Porter-Thomas distribution of r n /r n by a 
group of 10 equally probable r n 's per spin state. 

The value of r is assumed to be constant at 125 meV. It might be mentioned that 
125 meV lies within the probable error range of 57 of the 63 resonances measured by the 
Saclay group between 0 and 1000 eV. This assumption, that of constant r , is consis- 
tent with theory which indicates that should be distributed as x^ with a large num- 
ber of degrees of freedom (ref. 18). Since the mean square deviation from the mean, 

2 

for a normalized x distribution, is two over the number of degrees of freedom, r 
should remain relatively constant from resonance to resonance (ref. 19). 

A Doppler -broadening Breit-Wigner code is used to compute 10 or 20 resonance 
shapes, according to whether a merged spin-state or a separated spin-state cross section 
set is to be generated. The Doppler -broadening code computes values for the resonance 
cross sections from (E Q - 100) to (E Q + 100) eV by 0. 5-eV increments, providing 401 
energy dependent values of the cross sections per resonance. Figures 2 and 3 illustrate 
representative values of scattering and capture cross sections after Doppler -broadening 
for a merged spin-state calculation. The Wigner two-species level spacing distribution 
is chosen for the determination of the mock cross section's level spacings (ref. 17). The 
Wigner two -group distribution of level spacings shows the best agreement with the 
experimental data. 

A single-group Wigner distribution is assumed for the merged spin calculations. 

This assumption is used to evaluate the error introduced into merged spin-state calcula- 
tions. No correlation is introduced into the spin -dependent cross sections. Again, 
stratified sampling is used. The cumulative distribution function of the Wigner level 
spacing distribution is divided into 10 equally probable strata, and one shown in table HI. 
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s-Wave scattering cross section, 


Energy relative to resonance energy, |E 0 - E|, eV 

Figure 2. - Group s-wave scattering cross sections for gold. 
Potential cross section, 10 barns. 

Energy relative to resonance energy, |Eq -E|, eV 

Figure 3. - Group s-wave capture cross sections for gold. 
Potential cross section, 10 barns. 

TABLE m. - GROUP REPRESENTATION OF s-WAVE STATISTICAL 

DISTRIBUTIONS OF NEUTRON WIDTH AND LEVEL SPACING 


i 

24 keV FOR GOLD 



Population 

range 

Neutron width, 

r„(J = i), 
meV 

Neutron width, 

r n (J = 2), 

meV 

Level spacing, 
Dj(J = 1), 
eV 

Level spacing, 
Dj(J = 2), 
eV 



0 to 1. 0 

930 

930 

44.0 

24.0 



0 to 0. 1 

9.3 

9.3 

11.5 

6.5 



0. 1 to 0. 2 

40.5 

40.5 

20.0 

11. 0 



0. 2 to 0. 3 

102.3 

102.3 

26.0 

14.0 



0. 3 to 0. 4 

226.5 

226.5 

31. 0 

17. 0 



0. 4 to 0. 5 

339. 0 

339. 0 

38.0 

20.5 


| 

0. 5 to 0. 6 

539.9 

539.9 

43.0 

23.5 



0. 6 to 0. 7 

827.7 

827.7 

49.0 

27. 0 



0. 7 to 0. 8 

1246.2 

1246.2 

58.0 

31.5 



0. 8 to 0. 9 

1925. 1 

1925. 1 

68.0 

37.0 



0. 9 to 1. 0 

4043. 5 

4043.5 

95.5 

52.0 






These input data are useful for the stochastic generation of the spin -dependent cross 
sections. The procedure used to generate the cross sections in the unresolved energy 
region follows: 

(1) A random integer, uniformly distributed in the interval 1 to 10, is chosen. 

(2) One of the 10 resonances of a spin-state family, which corresponds to that random 
integer sampled in step 1, is chosen. 

(3) A random integer, uniformly distributed in the interval 1 to 10, is chosen. 

(4) Associated with the random integer just chosen is a level spacing, which is used 
to orient the next resonance obtained. 

(5) Steps 1 to 4 are repeated until the cross sections associated with 10 energy 
points have been generated. The same procedure is repeated for the other spin state 
with a new random number sequence to ensure stochastic independence. 

A question is posed as to how many of simulated cross sections generated, are 
necessary to represent adequately the entire set of one million cross sections. The 
question is answered in the following manner. The one million cross sections, for a 
single spin state, are divided into 2000 groups of 500 cross sections each, as well as, 
other various sized groupings. The mean value of each group is computed and the dis- 
tribution of these mean values are plotted. It was found that the distribution of the mean 
values for groups of 10 000 assumes normality. Thus, by the Central Limit Theorem 
(refs. 20 and 21) a group of 10 000 is an adequate sample size for stochastic independ- 
ence. The best set of 10 000 is chosen by matching the moments of the entire sample 
universe to each set of 10 000. The set whose first four moments best match the first 
four moments of the 10 cross sections is chosen as the most representative cross 
section set, for its particular spin state. This process is repeated for the other spin 
state. 

The procedure has generated sets of cross sections that are compatible with the low- 
energy neutron spectrographic data and the assumed statistical distributions associated 
with the levels and their spacings. A stochastically simulated total cross section set for 
gold is shown in figure 4. It is clear that, although the average total cross section is 
13. 67 barns, local values at 24 keV vary from 6 to 40 barns. 
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MONTE CARLO ANALYSIS OF GOLD SHELLS 


The Monte Carlo transmission code described and the stochastically generated mock 
cross sections together with parametric values of both the potential cross section and the 
p-wave capture cross section are used to calculate values of transmission for 24-keV 
neutrons. Values of the neutron transmission are computed for the various gold spheres 
of Schmitt and Cook (ref. 7) as well as for the gold sphere of Belanova et al. (ref. 8). 

The p-wave capture cross section is treated as an additive constant to the stochast- 
ically generated s-wave cross section sets because of the absence of measured p-wave 
statistics. Since there are four uncorrelated J-states for the p-wave resonances and 

o 

if the Porter-Thomas distribution for p-wave resonances applies, it is a \ distribution 

with two degrees of freedom. This distribution is not as widely a varying function as the 
o 

x with one degree of freedom (ref. 19), and does not emphasize the many small values 

ft 

of T n /r n as the x with one degree of freedom does. Thus the assumption of a con- 
stant p-wave capture cross section seems reasonable. 

With each value of potential cross section chosen, a new stochastic set of cross 
sections is generated and values of 6. 00, 8. 00, and 10. 00 barns are used as the potential 
cross section for gold. The potential cross section a is chosen as a parameter to 
conserve the independently measured value of the total cross section o^. for gold of 
13. 67 barns. For a given value of ct^, several values of the average p-wave capture 

cross section cr are introduced and values of the neutron transmission are computed 

P 

for each experimental shell. For each value of 0 p ot a value of <y Q is found which is 

compatible with the experimentally measured value of transmission. Figure 5 illustrates 



Figure 5. - Determination of p-wave capture cross section com- 
patible with experimentally measured transmission for Au I 
sphere. Potential cross section, 6 barns. 
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TABLE IV. - TRANSMISSION COMPATIBLE VALUES 
OF CROSS SECTIONS FOR Au I GOLD SPHERE a 


Potential, 

°pot 

s-Wave 

capture, 

% 

s-Wave 

scattering, 

% 

p-Wave 

capture, 

% 

Total, 

a tot 


Cross section, b 



6. 00 

0.499 

5. 151 

0.276 

11.926 

8. 00 

.499 

5. 147 

.205 

13.851 

10.00 

.499 

5. 150 

. 141 

15.790 


a Data from Schmitt and Cook (ref. 7). 
^See fig. 5. 


how the values of a are obtained for the Au I sphere of Schmitt and Cook (ref. 7) for 

P 

a of 6. 00 barns. This procedure is repeated for several values of a Table IV 

indicates the calculated values of cross sections that are compatible with the experimental 
transmission for the Au I sphere for assumed values of the potential cross section. How- 
ever, it is noted that the average total cross section ct^T in the last column varies. 

Since o^ 0 t is measured independently, the results in table IV can be interpolated for best 
values of a |., which is the only free cross section. In figure 6, a plot of against 

is shown. For this sphere transmission experiment, the potential and average 



Figure 6. - Determination of potential cross section compat- 
ible with experimentally measured value of total cross 
section. 
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p-wave capture cross sections that best match the measured value of the total cross 
section have been obtained as or j. = 7.80 barns and cF - = 0. 226 barn, respectively. 

P 

RESULTS AND CONCLUSIONS 

The Monte Carlo cross section generation process permits accurate compilation of 
the statistical aspects of neutron cross sections and their associated distribution func- 
tions. The Monte Carlo transmission code permits analysis of average values of trans- 
mission by neutron trajectories in the spherical shell experiments. 


Average Capture Cross Section 

A Monte Carlo analysis for the merged spin statistics gives significantly different 
results than a similar analysis using separated spin states. The results of such analyses 
are shown in table V for the Schmitt and Cook (ref. 7) spheres. For the merged results. 


TABLE V. - RESULTS OF GOLD SPHERE TRANSMISSION 
ANALYSES BY MONTE CARLO 
[ a Total cross section for all cases, 13. 7±0. 3 b.] 


Type of 
analysis 


Cross section, b 


s-Wave capture, 
% 

s-Wave scattering, 
% 

p-Wave capture, 
% 

Potential, 

°pot 

Unmerged 13 

0. 499±0. 012 

5. 15±0. 12 

0, 226±0. 025 

7. 8±0. 3 

Merged 13 

. 465±0. 010 

4. 45±0. 10 

. 194±0. 026 

8. 5±0. 3 

Merged 0 

. 485±0. 050 

3. 18±0. 40 

. 150±0. 015 

9. 9±0. 4 


a Value obtained from Duke University. 
^Strength functions from ref. 11. 
c Strength function from ref . 16. 


both the Columbia unseparated statistics and the merged Saclay data, are used to obtain 
average cross sections. For gold there is an increase of about 35 millibarns of unmerged 
over the merged results for the s-wave capture. This may be expected because the 
merged spin state uses a single -species Wigner distribution and forces level repulsion. 
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The separated spin-state analysis on the other hand, forces level repulsion only within a 
particular spin state, allowing the uncorrelated separate spin states to aline without 
introducing a spin-spin level repulsion. The average s-wave scattering cross section is 
likewise increased for the same reason. The value of the potential cross section is lower 
for unmerged spins due to the increased value of the average s-wave scattering cross 
section. The reason for the larger value of the average p-wave capture cross section, 
for the separate spin -state analysis, appears to be that, since the value of the potential 
cross section is lower, more capture is necessary to maintain the correct value of the 
transmission coefficient. 


TABLE VI. - AVERAGE CROSS SECTIONS FOR GOLD AT 24 keV 


Experiment 

Type of 

Shell dimen- 


Monte Carlo results 


Transmission, 


analysis 

sions. 






b 



cm 


s-Wave 

s -Wave 

p -Wave 

Capture 

Potential 






scattering 

capture 

capture 

cross 

cross 




Outside 

Inside 

cross 

cross 

cross 

section, 

section, 




radius 

radius 

section, 

section, 

section, 

V 

a pot’ 







°c~> 

^c~~’ 

b 

b 






s s 

c s 

C P 








b 

b 

b 




a Au I 

Merged 

7.62 

5.93 

4. 45±0. 10 

0. 465±0. 010 

0. 185±0. 025 

0. 650±0. 030 

8. 56±0. 31 

0. 876±0. 005 


Unmerged 

7.62 

5.93 

5. 15±0. 12 

. 499±0. 012 

.21U0. 026 

. 710±0. 031 

7. 80±0. 28 

. 876±0. 005 

a Au n 

Merged 

7.62 

5.08 

4. 45±0. 10 

0. 465±0. 010 

0. 204±0. 028 

0. 669±0. 032 

8. 54±0. 30 

0. 800±0. 004 


Unmerged 

7.62 

5.08 

5. 14±0. 12 

. 49 9 ±0. 012 

. 238±0. 029 

. 737±0. 033 

7. 78±0. 28 

. 800±0. 004 

b Long 

Merged 

3.55 

2.05 

4. 45±0. 10 

0. 465±0.010 

0. 259±0. 023 

0. 724±0. 029 

8. 49±0. 30 

0. 895±0. 002 

counter 

Unmerged 

3.55 

2.05 

5. 15±0. 12 

. 499±0. 012 

.275±0. 025 

.774±0. 031 

7. 74±0. 28 

. 985±0. 002 


^ef. 7. 
b Ref. 8. 


The final results are shown in table VI which compares the results of a merged 
Saclay statistics Monte Carlo analysis for the individual spheres with results for the 
separated Saclay statistics giving the individual cross sections derived from each sphere 
transmission experiment. The average unmerged value of for the Schmitt and Cook 
(ref. 7) experiments is 725±30 millibarns. The corresponding merged value of is 
660±30 millibarns. The average unmerged value of for the Schmitt and Cook experi- 
ments is likewise in good agreement with the recent gold activation measurement of Harris 
et al. (ref. 5) of 713±40 millibarns. The average unmerged value of a~ for Belanova 
et al. (ref. 8) is 774±31 millibarns and likewise the merged value is 724±29 millibarns. 


Average Potential Cross Section 

The results of the Monte Carlo analysis for gold indicate a potential cross section of 
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7. 80±0. 28 barns. This disagrees with a measured low-energy value of <Xp ot (11. 2 b, 
ref. 22) obtained by subtracting resonance contributions to the total cross section in the 
resolved region. The present result of 7. 80 barns at 24 keV may be converted to its 
low-energy value by a partial wave analysis using the rigid-sphere model (ref. 22) to 
provide a value of 8. 00 barns. With this model, the difference between CT pot at 24 keV 
and at low energy is small (2. 5 percent). Since <t , is given by 4 jt(R') 2 the radius 

POI _ ^ o 

implied by the potential cross section of 7. 80 barns is 7. 88x10 centimeter. This 

-13 

result is in remarkable agreement with the value of 7. 86x10 centimeter obtained from 
the strong interaction model R* = R = 1. 35xA^^ centimeter, where R' is the potential 
scattering radius and R is the nuclear radius. The discrepancy between the measured 
(ref. 22) a po ^ = 11. 2 barns and the smaller Monte Carlo result may be due to an under- 
estimation of the scattering cross section due to nearby low-energy resonances in ref- 
erence 22. The value of R' as given by the Monte Carlo analysis is compared with the 

197 

average radial density distribution for the Au nucleus as determined by electron scat- 
tering (ref. 24) in figure 7. The normalized density p(r)/p(0) is 0. 05 at a radius of 

-13 13 

7.80x10 centimeter and is virtually zero at 8. 00x10 centimeter. Therefore, there 

is good agreement between the Monte Carlo estimate for the nuclear radius and the nu- 
clear density distribution as ascertained by electron scattering experiments, whereas 
the value of R reported by Seth et al. (ref. 22) is considerably larger. 


R' region from Monte Carlo 
calculation (7.80±0. 20)xl0 -13 cm-v v 


R' region from Seth 
(9.40±0.10)xl(T 13 cm 



Radius, r, cm 


Figure 7. - Normalized radial density distribution compared with Monte Carlo potential scat- 
tering radius R 1 and R 1 of Seth et al. (ref. 23). 


The total average scattering cross section as determined by Monte Carlo is 
12. 95±0. 20 barns as compared with the measurement of Langsdorf et al. (see ref. 25) of 
12. 80±0. 20 barns; this indicates good agreement between the present Monte Carlo analy- 
sis and an independently measured value of a^. 
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It may be recalled that the only free parameters in the present analysis are the 

potential cross section and the average p-wave capture cross section. They are adjusted 

to fit the total cross section and the transmission, respectively. Thus the agreement of 

the potential cross section with the other independently measured and computed quantities 

and the agreement of the Monte Carlo total average scattering cross section with that of 

Langsdorf are additional contributory evidence for the validity of the Monte Carlo results. 

A method for the generation of cross sections in the unresolved energy region has 

197 

been illustrated and used to evaluate and interpret Au spherical shell neutron trans- 
mission experiments. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 20, 1969, 

129-02-04-04-22. 
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APPENDIX - METHOD OF STRATIFIED SAMPLING 


The production of the generated cross sections in the unresolved region entails the 
sampling from two probability distributions, the Wigner line spacing distribution and the 
Porter-Thomas distribution of neutron widths. Since both distributions are known and 
may be expressed analytically, one might be tempted to sample in a direct manner from 
both these distributions. One could continuously map the uniformly distributed random 
number generator on the two cumulative probability distributions and for each sampling 
of a random number in the interval (0, 1) obtain through the mapping, a value of the 
sampled variable (ref. 26). However, this process has been rejected for reasons which 
shall now be illustrated. The Wigner line spacing distribution may be written as 



If one makes the simple transformation, x = D/ (D) one obtains 


Pjt)=I -J^exp t-lx 2 
2 <D> V 4 , 


If one then forms the functional mapping 



and solves for X, one finds 

X = -LJ_ln(l - R ) 

v? 1 

Thus, to choose from the Wigner level spacing distribution in a continuous manner, 
first a random number must be chosen in the interval (0, 1), then that number subtracted 
from 1. Next the natural logarithm of the difference is taken, the sign of the natural 
logarithm is changed, the square root of the logarithm is taken, and that result is then 
multiplied by 2/y/W . Thus in the production of one random sample from the Wigner dis- 
tribution at least one subtraction, one multiplication, and three function evaluations are 
performed. Hence, for an IBM 7094-11, the time necessary for one sampling is of the 
order of 450 microseconds, a similar calculation for the Porter-Thomas distribution 
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gives much larger generation times. In order to reduce the amount of time required to 
generate cross sections and to reduce the variance associated with such a generated set, 
the technique of stratified sampling is used. 

Stratified sampling consists of dividing the population to be sampled into various 
strata with appropriate weighting and then the choosing values from each stratum to 
represent the distribution within that stratum (ref. 27). In this particular application, 
the cumulative probability distribution functions of both the Wigner and Porter -Thomas 
distribution are divided into 10 equally probable strata. Then the average value within 
each stratum is obtained. These values are chosen to preserve the average value of the 
entire distribution. The values, thus obtained, are checked using the Kolmogorov - 
Smirnov test for goodness of fit (ref. 28) and other statistical methods. Both the Porter - 
Thomas and the Wigner samples are accepted at the 99 -percent level of the Kolmogorov - 
Smirnov test. A comparison of the moments of stratified samples and some random 
samples is shown in table VII. This comparison indicates the power of the stratified 


TABLE Vn. - WIGNER LEVEL SPACING 


Sample 


Moments about mean 


size 

Mi 

** 2 

^3 

m 4 

^2 



Stratified sample 



10 

0 

0. 3209 

0. 1420 

0. 3403 

0.3265 

CO 

0 

.2732 

.0901 

.2423 

.2341 



Random sample 



40 

0. 0486 

0.3516 

0. 2266 

0. 5030 

0.7557 

100 

.0189 

.2972 

.1154 

.2873 

.3207 

500 

.0462 

.2846 

.1409 

.3053 

.3981 

1000 

0. 

.2732 

.0901 

.2422 

.2341 


sampling technique, in that the stratified sample size of 10 gives an approximation of the 
exact moments which is matched only after sampling randomly hundreds of values. The 
same is true for the Porter-Thomas distribution. 
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